Structure and heat capacity of Ne and Xe adsorbed on a bundle of carbon nanotubes 
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The structural and thermal properties of Ne and Xe gases adsorbed on the outer surface of a large 
nanotube bundle have been evaluated with computer simulation. The potential energy model and 
numerical techniques were used previously to study Ar [N. M. Urban, S. M. Gatica, M. W. Cole, 
and J. L. Riccardo, "Correlation functions and thermal properties of Ar adsorbed on the external 
surface of a bundle of carbon nanotubes", Phys. Rev. B 71, 245410 (2005)]. Heat capacity results 
for Ne and Xe exhibit peaks associated with reordering and "stripe" melting transitions for these 
gases. 

PACS numbers: 61.46.Fg, 68.43.De, 68.43.Fg 

I. INTRODUCTION 

In recent years there has been considerable discussion of the properties of gases adsorbed on bundles of carbon 
nanotubes. Measurements of thermal properties (isotherms, isosteric heats, and thermal desorption behavior) have 
provided information about the energies of these gases, while scattering and spectral data have provided information 
about their structures and dynamical properties ^&2££&iZJL In parallel with these experiments have been various 
theoretical and simulation studies, predicting virtually every property that is amenable to such measurements. Varying 
assumptions underlying the calculations pertain to the geometries and the adsorption potentials. 

Our group has undertaken a series of computer simulation studies of simple classical gases adsorbed on the outside 
of large bundles^ ' 10 ' 11 Since these bundle are large, the calculations actually assumed the bundle surface to be flat, 
which would be precisely correct if the bundle were actually infinite in radius or else if the bundle consists of a layer 
of nanotubes deposited on a surface. The potential experienced by an adsorbed molecule was evaluated by summing 
contributions from a line of parallel tubes, each of radius R; this approach neglects contributions from the interior of 
the bundle, an assumption that is plausible only for the first or second adlayer. Fig. Q] shows the potential energy that 
results for a Ne atom in such an environment, based on the potential model used in our group. That model ignores 
the chirality of the nanotubes by smearing out the carbon atoms into a cylinder made of continuous carbon matter. 
These assumptions are semiquantitatively accurate, at best. 

Our simulations employ either canonical or grand canonical Monte Carlo methods. In early work, the familiar 
Metropolis algorithm was used, but our most recent studyii of Ar (denoted I) employed the efficient Wang-Landau 
algorithm 12 , discussed below. The adsorption behavior derived from the model consists of a sequence of "striped" 
phases, with the number of stripes increasing as the gas pressure P is increased. The first of these, appearing at low 
P because of the very attractive potential, is the so-called "groove" phase, in which the gas is confined to a quasi-one 
dimensional vicinity of the potential minimum, nestled between two adjacent tubes. As P increases, a transition to 
a "three stripe" phase occurs, followed by the appearance of a well defined monolayer. An intriguing "second layer 
groove" phase has been seen in both simulations and experiments. Multilayer films can occur, but these have not 
been investigated because desorption precludes their observation, except at very low temperature. 

This paper focuses on the thermal properties of Ne and Xe, extending the Ar study, I. The next section summarizes 
the methodology and reports results for a very low density phase, called N = 1. The number N is the number of 
adparticles confined within one unit cell of the simulation (which is periodically replicated across the surface) . The 
extent of this cell in the lateral (x) direction is 1.7 nm, taken to be the spacing between (10,10) tubes. In the y 
direction (parallel to the tubes' axes), the extent is lOcr, where a is the Lennard- Jones diameter of the particles; the 
other gas parameter entering the calculation is e, the well depth of the adatom-carbon pair interaction. In the z 
direction, perpendicular to the bundle surface, the simulation cell extends a distance H = 4 nm, where a hard wall 
terminates the simulation volume. Since this boundary underestimates the (macroscopic) experimental gas volume, 
its presence tends to suppress desorption relative to what is seen in experiments. 

Sections Mil and IIVI describe simulation results obtained for Ne and Xe, respectively. Section [V] summarizes our 
principal results and draws conclusions. 

II. CALCULATIONS FOR N = 1 

The Wang-Landau algorithm is a simulation technique that generates the density of states g(U), the number of 
states per unit energy at energy U, for a specified number (N) of particles; we omit a descriptor N from g(U) for 
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FIG. 1: Contours of constant potential energy for a Ne atom (a) above the nanotube surface compared with (b) energy contours 
generated by the simple model described in the text. Shown are isopotential lines for V = —700, —600, —500, —400, —350, 
—300, —250, and —200 K, starting from the groove. The rapid change in area between curves for —400 K and —350 K (dashed 
lines) manifests a high density of states, giving rise to a heat capacity peak. 



notational simplicity. The procedure used to compute g(U) is discussed in I and in greater detail in a recent thesis 
reporting this workji 3 - From this function, thermal properties may be computed at any temperature T — l/(fce/3). 
Since classical systems have a specified kinetic energy per particle, 3/(2/3), we focus on the potential energy and 
omit that contribution except when showing numerical results for the heat capacity C(T). For example, the (mean) 
potential energy (U(N,T)) of the system is 

(U(N,T))= J Ug(U)e-? u du/[J g(U)e-? u du]. (1) 

The specific heat (in Boltzmann units) is then computed from 

C/(Nk B ) = 3/2 + [dU/d{Nk B T)] N . (2) 

Results for g(U) for N = 1 are shown in Fig. [^i 3 - This function exhibits a striking variation with U , which can be 
understood with simplifying models, as follows. For classical systems, g(U) dll is strictly proportional to the volume of 
phase space confined between surfaces of energy U and U +dU. This quantity, for N=l particle, was called the "volume 
density of states" in I, where it was evaluated from histograms of the potential energy distribution in real space; the 
result agrees with that found directly with the Wang-Landau simulation method. The lowest energy "states" originate 
from the vicinity of the groove region. The resulting low U behavior, g(U) ~ constant, was derived in I by assuming 
that the constant energy contours are circles, centered on the groove; ellipses yield the same analytical form. The two 
most striking features of g(U) are a sharp peak, near a potential energy U = V m = —363.373 K of and a divergence 
near the desorption threshold, {7 = 0. The latter divergence is of the form g(U) ~ |J7| -4 / 3 , where the power 4/3 was 
derived in I from the asymptotic van der Waals interaction, V(z) ~ — 1/z 3 . The divergent behavior just above the 
peak at V m , g(U) ~ (U — Vm)^ 1 ^ 2 was derived in I, by analyzing an assumed potential energy near the cylindrical 
surface at which the monolayer forms: V = V m + (k/2)(z — z m ) 2 . Here, the value V m is the holding potential of a 
particle near a single cylindrical tube; the quadratic functional form is analogous to that found near the equilibrium 
position on a flat surface. To understand the divergent behavior of g(U) just below the peak (U < V m ), we have 
invented a simple model potential that roughly characterizes the spatial dependence near the groove and monolayer: 

V mode i(x,z) =V m - Aexp(-a; 2 ) + (fc/2)(z-z m ) 2 . (3) 

Here A = V m and k — 200 K/A 2 are fitting parameters determined from the potential. This form is shown in Fig. [1] 
to resemble that of the "true" potential near the divergence for U ~ V m . Even for this simple model, we cannot 
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FIG. 2: Density of states g(U) for N = 1 Ne atom (full curve). This is compared with the prediction (dashed curve) derived 
from the simple model, Eq. [3] 

determine g(U) analytically, unfortunately, but we can evaluate it numerically. For the "true" potential the density 
of states was calculated directly using Wang-Landau Monte Carlo, and for the simple model it was calculated by 
discretizing space and binning the volume elements with potentials between U and U + dll to produce a histogram. 
As seen in Fig. [21 there is a qualitative agreement between the two curves for g(U) in the region of the divergence, 
just below U = V m . The divergence originates because there is a large volume, near the outermost ellipse surrounding 
the groove, for which the potential energy is close to this value. 

The specific heat of Ne and Xe, for N = 1, is analogous to that found for Ar. At low T, the value 5/2 is expected 
for the specific heat C j (Nks) because particles have two degrees of freedom due to potential energy of the confining 
potential within the groove, plus three more degrees of freedom from the kinetic energy. The peak in C(T) near 
T = 80 K for Ne, or T = 280 K for Xe, is due to excitation out of the groove, on to the cylindrical surface of 
the nanotube. This is confirmed by model calculations, as described in I, in which particles are located either in 
the groove region or on a nominally flat surface. The high T limit, C/(Nks) = 2, is attributable to the potential 
energy contribution from one degree of freedom, associated with motion perpendicular to the surface. This behavior 
is depicted, along with higher coverage results, in Figs. [3] and 01 



The Wang-Landau method was used to calculate the heat capacity C(T) for higher coverages (Fig. [3]).— For Ne, 
the groove phase occurs at a density of N = 9 adatoms in the simulation cell; the three stripe, monolayer, and second 
groove phases occur at TV = 27, 72, and 81, respectively. 

The low density N — 1 results show a heat capacity peak near T — 80 K, due to excitation of Ne atoms out of the 
groove. The higher coverage results have a similar peak (referred to here as the second peak), as well as an additional 
low temperature peak (first peak), i.e., a characteristic double peak structure in C(T). 

Monte Carlo configurational samples indicate that the second peak is due to substantial excitation of adatoms, as 
in the N — 1 case. The excitation is primarily from the groove into the stripe or monolayer regions. The second peak 
occurs at somewhat lower temperature than in the N = 1 case, because it is easier to excite atoms out of a densely 
packed groove than a nearly empty one. This is a result of repulsive mutual interactions that predominate when the 
groove nears completion. 

The lower temperature, first peak in the heat capacity, in contrast, is not associated with promotion of adatoms 
into different regions of the substrate. Rather, it is due to "energy broadening" , in which the adatoms remain at or 
near the same sites, but are thermally excited over a broader range of energies. This phenomenon is exhibited in 



III. NE ADSORPTION AT FINITE COVERAGE 
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FIG. 3: Specific heat of Ne for various coverages N (explained in the text), as a function of T. 
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FIG. 4: Ne density distribution in the three stripe (N — 27) phase at various temperatures (cross sectional view of adjacent 
nanotubes). The T — 20 and 75 K panels show the Ne distribution at the first and second peaks in the heat capacity. 



Fig. |4]for the three stripe phase. The distribution of Ne atoms indicates that at the first peak (T = 20 K), most of 
the atoms remain in the stripes, albeit with some promotion into the full monolayer. The atoms remaining in the 
stripes spread out, confined less strongly in the stripe potential minima, as a result of energy broadening. At the 
second peak, atoms in the stripes have been fully promoted into the monolayer and the stripes disappear. 
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FIG. 5: Specific heat of Xe for various coverages N, as a function of T. 



IV. XE ADSORPTION AT FINITE COVERAGE 



The Xe behavior, shown in Figs. [5] and [6l is qualitatively similar to that of Ne, as might be expected. There are 
some differences of detail, however. For example, the peaks are smaller and occur at much higher temperature in the 
Xe case. The higher temperature is attributable to the higher energy scale associated with the interactions' strengths; 
the gas-gas well depth for Xe is a factor of seven higher than that of Ne and the adsorption potential is nearly a 
factor of three more attractive. (The behavior of C{T) of Ar, explored in I, is intermediate between these two.) We 
surmise that the smaller peak height is due to the fact that the importance of the adhesive interaction, relative to 
the cohesive interaction, is reduced in the Xe case, so the peak in C(T), attributed primarily to the interaction with 
the substrate, is less narrow for Xe. For the same reason, we believe, "stronger" cohesive interactions, arising from 
subtle intra-film interactions, produce wiggles for C(T) of Xe, not found for Ne. (We have confirmed that these are 
not artifacts of the calculation by refining our simulations.) We conjecture that the intermediate "wiggle" in C(T) 
for the monolayer film, near T = 150 K, is due to the two non-groove stripes melting into the monolayer (compare 
the T = 100 and 150 K density distributions in Fig. [5]); the lower-energy groove stripe itself remains intact until the 
groove-to-monolayer promotion at the higher temperature peak near T — 250 K. 



V. CONCLUSIONS 



We have derived the evolution of density with T and the heat capacity for Ne and Xe gases adsorbed on an array 
of nanotubes. The behavior is qualitatively consistent with that reported previously in I for the case of Ar. These 
predictions, along with those of the isotherms ; 9 ' 10 ^ 1 are testable experimentally and we look forward to such results, 
from which conclusions may be drawn. The principal assumptions that may or may not be validated by such studies 
pertain to the geometry and the interactions. The latter type of assumption would seem more reliable, since the 
interactions are semi-empirical, derived from experiments involving adsorption on graphite. However, extrapolation 
to the nanotube case, as done here, represents a further, simplified assumption: that the chemistry of the atoms is not 
sensitive to their structure. While this might seem plausible for C, it remains to be verified. More severe, presumably, 
is the assumption that the interaction arises from a smeared-out array of C atoms. While clearly open to question, 
this assumption might be justified by the fact that since the assumed pair potential is already a rough approximation, 
the continuum approximation is not much more drastic. Of course, possible commensurate structures are not found 
in the continuum approach. 

Concerning the geometry, the assumption of a unique radius for all tubes is another oversimplification^ as is 
the presumption that the bundle radius is essentially infinite. While more detailed studies are feasible, they require 
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FIG. 6: Xe density distribution in the monolayer (N = 45) phase at various temperatures, shown as cross sectional view of 
adjacent nanotubes. The T = 100 and 250 K panels show the Xe distribution near the first and second peaks in the heat 
capacity. 



significantly greater computational effort. 
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